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Abstract 

We write down matrix models for Ising spins with zero external field on the vertices of dynamical 
triangulated random surfaces (DTRS) and dynamically quadrangulated random surfaces (DQRS) 
and compare these with the standard matrix model approach which places the spins on the dual (p^ 
and 0* graphs. We show that the critical temperatures calculated in the DTRS and DQRS models 
agree with those deduced from duality arguments in the standard approach. Using the DQRS model 
we observe that the Ising antiferromagnet still undergoes a phase transition to a Neel (checkerboard) 
ordered ground state which is absent because of frustration in the other cases. 
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The work of Boulatov and Kazakov [y solving the Ising model on dynamical surfaces predated both 
the recent explosion of interest in matrix models and the use of conformal field theory methods for 
c < 1 models coupled to 2d gravity They used techniques developed in earlier work by Brezin et.al. 
and Mehta Q to cast the Ising model on dynamical 0^ and (j)^ graphs (ie the duals of DQRS and DTRS 
respectively) as two matrix integrals. The partition function of the Ising model on some set of random 
graphs G" with n vertices can be written as 

z„(/3, if) = 51 E E ^S'^^^j + ^ E (1) 

G" {a} \<y> « / 

where is the connectivity matrix for a given graph. Summing over the number of vertices gives 



Z(c,.9,i/) = ^^(^^^-|^j Z„(/3,i/) (2) 

where c = exp(— 2/3). Boulatov and Kazakov showed that Z was given by the free energy of two N x N 
hermitian matrices U and V which in the case of 0^ graphs with external field H = was: 

F^.{c,g,H) = ^log (^J d^''M,d^'M,e^p-trj (^M + M _ ,m,M, + -^Mf + ^M^^^ (3) 
and for (j)^ graphs: 

F,.{c,g,H) = ijlog i^J d^'Ahd^'M,e^p~tr^ + ^ " ^^^i^^^^ + l^^' + l^^^)) ' (4) 

In these equations the Mi matrices can be thought of as representing <t = +1 vertices and the M2 matrices 
(7^—1 vertices. Rescaling to obtain (j)"^ and (j)^ couplings of g/N and g/VN respectively shows that the 
correct topological factor of ^ — E' + F is obtained for each graph, where V is the number of vertices E 
is the number of edges (propagators) and F is the number of faces. Both the 0* and (p'^ models displayed 

a third order phase transition with critical c values of c^^^j = 1/4 and c^^-j = (\/28 — l)/27. The critical 
exponents were the same for both lattices as one might have expected from universality and satisfied the 
hyperscaling relations, but were different from the fixed lattice Onsager exponents. The values of these 
new exponents for the dynamical case were confirmed by the Liouville theory calculations of KPZ and 
DDK in § in the light -cone and conformal gauges respectively. 

The same universality that gives identical exponents on (fi^ and (j>'^ lattices might be expected to 
carry over to DTRS and DQRS |^ and numerical simulations Q provide some confirmation of this, giving 
satisfactory agreement with the predictions of ||^, |^ for both DTRS and (jy^ graphs. It would be reassuring, 
however, to write down directly matrix models for Ising spins on DTRS and DQRS instead of relying 
on duality arguments to find the critical temperature and universality to find the exponents. We have 
another motivation for considering such models: An Ising anij-ferromagnet will display a phase transition 
to a Neel ordered ground state on "loosely packed" lattices where only even loops are present This 
is clear from the strong coupling expansion of equ.([^) 

Zv{(3, 0)=Y1 2^(cosh(/3))^ Yl (tanh(/3))'^"^'*^ (5) 

loops 

where the sum runs over closed loops with each edge traversed only once. Equ.(|^) possesses a /3 — *■ 
— /? symmetry if all possible loops are even in length, which translates into a mapping between the 
ferromagnetic (/3 positive) and antiferromagnetic (/? negative) partition functions and hence a phase 
transition in the antiferromagnet. The and (j)^ Ising models will contain loops of both even and 
odd length so frustration will suppress the antiferromagnetic transition. The f3 — > —(3 symmetry which 
corresponds to c ^ 1/c should hold at the level of the partition function or free energy and it is clear from 
equ.(|) and equ.(^ that the 0* and (f)^ matrix models do not possess this symmetry. An antiferromagnetic 



-"^This is not entirely trivial as it appears that situating spins in the centre of simplices in the 3d equivalent gives very 
different results from situating spins on the vertices ]7l 
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transition will also be absent on DTRS, but for DQRS we would expect to see a transition as all loops 
must be of even length. The matrix model for Ising spins on DQRS should therefore possess a /3 — > —(3 
symmetry. 

To write down the appropriate matrix models we change our perspective slightly and think of the 
potential terms in the matrix integral as representing the faces of the graph rather than the vertices |^ . 
For a one-matrix integral with action of the form tr[(j)^ /2 + {g / y/N)(j)'^) this still gives the correct counting 
as each internal vertex gives a factor of N and each face (triangle in this case) gives a factor of TV- ^ 
give a total of N^^^/"^ . As iF = 2E this produces the correct factor of 0. Similarly, an action of the 
form tr{4>^ /2 + {g/N)(j)'^) which generates squares will give N^~^ and in this case 4F — 2E, so we again 
find . Let us now apply similar considerations to the Ising model on DTRS with the external field H 
set to zero for simplicity. The triangular faces will be composed of two types of edges, S with the spins 
at both ends the same, and D with the spins at both ends different. Each S edge will contribute a factor 

cosh(/3) + sinh(/3) = cosh(/3)(l + c*) (6) 

and each D edge will contribute a factor 

cosh(/3) - sinh(/3) = cosh(^)(l - c*) (7) 

where we have introduced c* (1 — c)/(l + c), which is the transformation effected on c by taking 
the dual transformation of the temperature f3* ~ — l/21ogtanh(/?). There will be two possible types of 
triangles in the model, one with all the spins the same giving a trS^ term in the action, and one with 
two spins the same and one different giving an trSD'^ term. We choose to put the factors in equs.(^,0) 
in the propagators so the matrix model action we arrive at is 

UoTRS^-trf- —L^_^S'+ d'+s'/3 + SdA (8) 

g \2cosh(p)(l + c*) 2 cosh(/3)(l — c*) / 

where S and D are again N x N hermitian matrices. 

We can show that this is an Ising model with a critical temperature agreeing with that deduced from 
duality by transforming it to an 0(1) representation of the Ising model derived by Eynard and Zinn- Justin 
§. A rescaling S 5/cosh(/3), D ^2(1 - c*)!)/ cosh(/3), gcosh^{/3) g gives 

Udtrs = -tr (d\1 + 2S{\ - c*)) + \ ,, S'' + S^/s) . (9) 
9 \ 2(1 + c*) J 

The first term can be simplified at the expense of introducing a more complicated set of interactions 
by defining Si ^ 1 + 25(1 - c*) which gives, on scaling Si/{1 + c*), 2^(1 - c*f{l + c*)^ g, 

2D^{1 — c*)'^(l + c*)2 ^ and dropping constant terms 

t,,,,,.^t.(i,^a + g-fs;+ '^''-^]" + ''' 5.). (10) 

We have carried out these manipulations because it was shown in [|| that defining A — (Mi — Af2)/2, 
S = M2 + M2 + 1 + c transformed the (j)'^ action in equ.(|^) into 

which is identical in form to equ.(^). Performing a gaussian integration over A in this action allowed 
saddle point methods to be used to determine the critical behaviour of the model in the large N limit. 
It was shown that the resolvent of the matrix S 

Mz) = i,tr-^ (12) 



was given by 



N z- 



u;o{z) = j^i^V'iz) - V'{-z)) + ^c.(z) (13) 



^1 would like to thank Jan Ambj0rn for sorting out my confusion on this point 
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where we have written the action as A^S + V{S) and uj{z) — e a a;+(z) + e^" uj-{z) with, for a cubic 
potential V{z) 

1 /3 

u;±{z) ^ (^^/l~^¥/^Tib/z'j (^-/T^WJT^ ± iba / z^j . (14) 

For a potential of the form W{S) ~ S'^ /"i — aS'^ + f3S the parameters in the solution b and a were related 
to the parameters in the potential a and /3 by 

6(3f7-l) 6\/3a 

62(3(7-5) = 9/3. (15) 
The critical Ising model corresponded to <t = 1 which gave, on eliminating b, 

6a^+P^ = 0. (16) 
If we substitute the values for the cj)'^ model {a — 2c, (3 = (3c— l)(l+c)) into equ. ([l^) we get 27c2 + 2c— 1 — 

,3 

0, which reproduces the correct value of c^.^^^ = (v28 — l)/27. For the DTRS model we replace c by c* 
which on substituting into equ.(|l|) gives 21{c*Y + 2c* - 1 = 0. The critical value of c* on DTRS is thus 
(V28 - l)/27 so the critical value of c (= (1 - c*)/(l + c*)) is c^f^^^ = (14 - v^)/(13 + ^/^), which is 
the dual of the value on (^•^ graphs. It is also worth remarking that /3 —(3 transforms c* — c* which 
is obviously not a symmetry of equ. (p^) , so there will be no antiferromagnetic transition on DTRS. 

To write down the Ising model on DQRS we use similar notation, representing same spin edges by 
5 and different spin edges by D. In this case we will attach the weights to the interaction terms that 
generate the various possible square faces, though it is also possible to attach the inverse factors to the 
propagators. The action is given by 

UoQRS = ^tr + + + + ^-^{SDSD + 2S'd')^ (17) 

where 

.91 = cosh'(/3)(l + c*)2 

52 = cosh2(/3)(l-c*)2 

51 - cosh2(/3)(l + c*)(l-c*). (18) 

These couplings are the square root of those one might naively write, but it should be remembered that 
each edge is shared by two squares. Using 

cosh(/3)(l ± c*) = J (19) 
cosh(/5)(l q= c*) 

and scaling S S cosh{0), D —> I?cosh(/3), g — ^ gcosh^{0) we can write this as 

= J'^ {I'' + l^' + 4(1^^' + 4(TT^^' + 2(l-c*)(l + c*) (^^^^ + 

(20) 

We can massage the cj)'^ model of equ.(|^) into a similar form by defining A = {AIi + Af2)/V2, B = 
(Ml - M2)/V2 and then rescaling A V2A/\/]~~c, B V2B/y/l + c, g ^ 2g to get 

[/.4 = —tr ( -A^ + -B^ + —^—A* + — J— -S* + -(ABAB + 2A^B^)] . (21) 

g \2 2 4(1 -c)2 4(l + c)2 2(l-c)(l + c)^ / 

which is identical to equ. ( pO| ) up to the replacement of c* by c. We know that in equ. ( pT| ) represents 
an Ising model with c„.jj = 1/4, so Udqrs in equ. (20) must represent an Ising model with c*^^^ — 1/4. 
This gives c^^^"^ = {X~ c*rit)/(l + c*rit) = 3/5, which is identical to the value deduced from duality and 
the (j)^ model result. In addition we can see that Udqrs possesses an obvious c* ^ — c* symmetry if one 
also exchanges S and D so we can deduce that the Ising antiferromagnet on DQRS will display a phase 
transition. One might worry that the summation over quadrangulations would disrupt the Neel order 
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of the low temperature phase, casting doubt on the above observation. We can see, however, that this 
is not so by considering the "flip" moves that implement the sum over quadrangulations in simulations. 
The two possible flips for squares maintain a connection between unlike spins and thus cost no energy in 
a Neel ground state - see Fig.l. 

To summarize, we have written down matrix models that represent Ising spins on the vertices of 
DTRS and DQRS and shown that their critical points agree with those deduced from duality for Ising 
spins on (p^ and (j)^ graphs. Ising spins on DQRS also manifest an antiferromagnetic transition because 
only even loops are possible in this case. It would be interesting to include a non-zero external field in 
the models in this paper. It is perhaps worth remarking in closing that complex matrix models generate 
surfaces composed of only even loops ]lO| so one might also consider using a two complex matrix model 
to examine the behaviour of the Ising antiferromagnet. 
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Figure Caption 



Fig. 1. The two possible flip moves on adjacent squares in DQRS both preserve the "checkerboard' 
order of the ground state. Up spins are shown as solid dots and down spins as open dots. 
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